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We experimentally analyzed the statistical errors in 
quantum-state estimation and examined whether their lower 
bound, which is derived from the Cramer-Rao inequality, can 
be truly attained or not. In the experiments, polarization 
states of bi-photons produced via spontaneous parametric 
down-conversion were estimated employing tomographic mea- 
surements. Using a new estimation strategy based on Akaike 's 
information criterion, we demonstrated that the errors actu- 
ally approach the lower bound, while they fail to approach it 
using the conventional estimation strategy. 

PACS numbers: 03.67.-a, 42.50.-p, 89.70,+c 

One of the central features of quantum mechanics is 
that it does not allow obtaining complete information 
about an individual quantum system without errors [jjj. 
The Holevo bound and the no- cloning theorem are promi- 
nent manifestations of the restrictions on acquiring in- 
formation from a quantum system Although it is 
possible to estimate all aspects of a quantum state by 
performing a series of distinct measurements on identi- 
cally prepared particles, that is, by tomographic mea- 
surements [HHH , there is still a lower bound on the sta- 
tistical errors. Therefore, it is important to establish 
an estimation strategy that can attain the lower bound 
to improve the accuracy and sensitivity of the precision 
measurements that exploit quantum states Hff]||] and to 
develop quantum information technology that requires us 
to faithfully prepare quantum states gj. 

In this article, we report our experimental analysis of 
the statistical errors in quantum-state estimation. Vari- 
ous polarization states of bi-photons produced via spon- 
taneous parametric down-conversion were estimated em- 
ploying tomographic measurements Although esti- 
mating such a two-qubit system has already been in- 
vestigated in detail by James, et al. ||, achieving the 
lower bound on statistical errors, which is derived from 
Cramer-Rao inequality, is highlighted in this article. We 
experimentally examined whether the lower bound was 
truly achieved, and found that there are situations in 
which the use of redundant parameters to estimate the 
quantum state leads to ambiguities in the estimates and 
prevents the errors from approaching their lower val- 
ues. We then developed a new strategy for estimating 



the quantum state based on Akaike 's information crite- 
rion (AIC) ||, which provides a rigorous way of selecting 
the most suitable number of parameters characterizing a 
quantum state. Using the new estimation strategy, we 
found that the errors in the quantum-state estimation 
actually approach the lower bound. 
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FIG. 1. Experimental setup for producing various polariza- 
tion states of bi-photons and measuring them. 

A rough sketch of the experimental setup is shown 
in Fig. |]. We generated bi-photons at a wavelength 
of 532 nm via type-I spontaneous parametric down- 
conversion (SPDC) with two beta-barium-borate (BBO) 
crystals attached together jl0| . The method used to pro- 
duce various polarization states of bi-photons was almost 
the same as that used in Refs. [| 11 ; the major differ- 
ence was that we used an ultra-short pulse laser (the third 
harmonics of a mode-locked TbSapphire laser) as a pump 
beam for the SPDC [0. We produced three quantum 
states, almost pure separable state (APSS), highly entan- 
gled state (HES), and very noisy mixed state (VNMS), by 
adjusting the pump-beam polarization with a half-wave 
plate by modifying the relative time delay between 
the horizontal and vertical components of the pump beam 
with a pre-compensator 12 1, and by inserting decoherers, 
i.e., two de-polarizers, into either paths of the down- 
converted photons ||ri|| . The coincident (within 6 ns) 
detection events, n v , on both single-photon detectors 
(HAMAMATSU H7421-40) at a certain polarizer's set- 
ting (i.e., projector) \m v )(rn v \ (which was determined by 
the half-wave plate, quarter-wave plate, and polarizer on 
each path of the converted photons) were counted by us- 
ing a time interval analyzer (YOKOGAWA TA-520) dur- 
ing data acquisition time t. To eliminate ambient pho- 
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tons, interference filters (FWHM: 8 nm) were used. The 
typical coincidence counting rate was about 500 Hz. 

The polarization states of the bi-photons (two-qubit 
states) can be characterized by a 4 x 4 density matrix with 
15 independent real parameters, pe = TqTq/Tt[T@Tq], 
where Tq is the complex lower triangular matrix with 16 
real parameters = {©} P^l- To determine 

these parameters, coincidence counting measurements 
at 16 particular polarizers settings, {|m„)(m„|}*=i [@i 
were used as the tomographic measurements p). From 
these coincidence counts, {rvlili = {N}, parameters 
{0} and thus quantum state pe were estimated using 
maximum likelihood estimation (MLE) |Jl3||l|,[l6| . 

We assessed the errors in the estimation as fol- 
lows. If the estimation procedure is repeated r times 
under identical conditions, we have r slightly differ- 
ent quantum states, {pQ.Yi=\^ due to statistical er- 
rors. Let the true quantum state be pe i where 
{0o} is the true values of the parameters. Then 
we calculate the average Bures distance J2l],|2^] be- 
tween the true state and the estimated states, i.e., 
2(1 - F(p eo ,Pe)) = 2(1 - £ ELi FjPfh, P^)), where 

/pe;]) is the fidelity §. 



F(pe ,P4.) (= Tr^^/peoPe., 
In the experiments, we repeated the estimation procedure 
nine times (i.e., r=9). The estimated quantum state us- 
ing all the data for the nine trials, i.e., IE,— 1 n v = 
{N }, is assumed to be the true state, pe [0- Keeping 
the experimental conditions identical, the average Bures 
distances were measured for five data acquisition times 
(t = 0.2 s, 0.5 s, 1 s, 2 s, 5 s). 

Experimentally measured average Bures distances (di- 
vided by 2) between the estimated states and their true 
states (solid points) and those obtained by Monte Carlo 
simulation (empty points) are shown in Fig. || (a) as a 
function of the data acquisition time multiplied by nui- 
sance parameter C, which corresponds to the coincidence 
counting rate without any polarizers (i.e., |m„)(m„| 
/ [Q. The simulations were carried out by artificially 
producing the 16 coincidence count data {N} accord- 
ing to the aforementioned true states. We estimated 
the states from the simulated data and repeated this 
procedure 200 times (i.e., r=200). Therefore, the val- 
ues presented in Fig. (empty points) correspond to 

1 — 2S0 Ej=i F(PQo 5 Pq ) ■ The results of the experiments 
and numerical simulations are in good agreement, even 
though the errors in experimentally estimating quantum 
state may stem not only from a statistical origin, but also 
from systematic errors in the experiment itself Jl9| . 

Next we examined whether the theoretically expected 
lower bound (20) on the Bures distance was attained or 
not in each experiment. The lower bounds are shown in 
the inset of Fig. || (a). They can be derived as follows jl7) . 
Suppose that the estimated state is in the neighborhood 
of the true state, so that the Bures distance between them 
can be written as [fn|,El2lM|2ll 
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where [■/y(O)] = J s (0) is the so-called symmetric loga- 
rithmic derivative (SLD) Fisher information matrix |2§| l 
and [V ij ] = [(0*(N) - 6 l )(6 j (N) - 6 j )] = V is the covari- 
ance matrix. This means that the Bures distance can 
be locally viewed as a distance on a Riemannian mani- 
fold equipped with a metric structure defined by the SLD 
Fisher information matrix [EmJUIlH ' 
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FIG. 2. Average Bures distances (divided by 2) be- 
tween true states and states estimated by MLE (a) and by 
MAICE (b), as a function of data acquisition time t multiplied 
by nuisance parameter C (« 500). The solid points (vertical 
error bars correspond to one standard deviation) represent 
experimental results and the empty points represent numer- 
ical results. The insets show the asymptotic lower bounds 
explained in the text. 

The Cramer- Rao inequality p6| , p7t provides an asymp- 
totic lower bound (CR bound) on covariance matrix V 
in Eq. ([[]): V > J _1 (0)le=eo, where J(9) is Fisher's 
information matrix p9| . The CR bound tells us that no 
estimation strategy exceeds this bound in the asymptotic 
region. With a little calculation, we have J(O) oc Ct [ p9[ . 
Therefore, we can define J(O) as J(6) = CtJ(Q). Conse- 
quently, from Eq. (Q) , the asymptotic lower bound on the 
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Bures distance between the true state and the estimated 
state is given by 

2(l-i?(p eo ,/>e))> i ^Tr[J s (e )(J- 1 (e)|e=eo)]- (2) 

We can see that the lower bound decreases as inversely 
proportional to the data acquisition time as shown in the 
inset of Fig. || (a). 

It is known that the CR bound can be asymptotically 
achieved by using MLE [|27|j30| ; however, the experimen- 
tal data in Fig. ||| (a), especially those for HES and APSS, 
show that the average Bures distances did not decrease 
as inversely proportional to t, so the CR bound was not 
achieved in these experiments. This means we need an 
estimation strategy other than MLE. 

To clarify the necessity of a new strategy, we briefly 
review the estimation strategy used in the above exper- 
iments. We define the true probability density function 
(PDF) obtaining the 16 coincidence counts {N} (with 
t = 1 s for each) as P{N\Q )\ t=1 = P (N)\ t=1 and 
its parametric model as P(N\Q)\ t= i. The relative en- 
tropy (|^J2^] between them is 

D(P (N)\\P(N\e))\ t=1 

oo oo 

= E- E (Po(N)HPo(N)) 

m=0 niQ— 

-P o (iV)ln[P(iV|0)])| t=1 , (3) 

and it takes a positive value, unless Po(N) = P(N\Q) 
for all {N} (in this case, D{P (N) \\ P(N\e))\ t= i = 0). 
Therefore, by minimizing the relative entropy in Eq. (j^) 
with respect to {O}, we can find the point {0(iV)} closest 
to the true point {Go} in the 16-dimensional parameter 
space. Because the approximation 

oo oo 1 

E • • • E W MP(iV|6)])| t=1 « - ln[P(JV|0)] (4) 

m— nie— 

is valid as the data acquisition time t is increased 
infinitely, this minimization is accomplished by max- 
imizing the logarithm of the likelihood function (log- 
LF) ln[P(7V|0)]. Namely, by finding parameters 
{0(A)} that maximize log-LF, the quantum state 
can be estimated (maximum likelihood estimation or 

MLE IJUJHII) S MLE has been used for estimat- 
ing the quantum states j|,[|[l^]. However, the ill-behavior 
of MLE in estimating the nearly degenerate states, such 
as HES and APSS, is clearly shown in Fig. || (a). 

The ill-behavior can be explained as follows. Quan- 
tum state pq has been implicitly assumed to be non- 
degenerate (i.e., in the interior of the space of the density 
matrix). Therefore, if the true state is a (nearly) degen- 
erate state, such as HES or APSS, (i.e., in (the neighbor- 
hood of) the boundary of the space of the density matrix, 
where one or more eigenvalues vanish), the density ma- 
trix can be completely characterized using less than 15 



parameters. Subsequently, the surplus parameters give 
rise to several local maximums of the log-LF in the 16- 
demensional parameter space and prevent the errors in 
the estimation from approaching their lower values. 

Thus, for attaining the lower bound, finding the ap- 
propriate number of parameters seems to be requisite. 
The Akaike's information criterion (AIC) || provides a 
rigorous method for eliminating redundant parameters. 
It is defined by 

AIC {k) (6) = -2 x ln[P (fe) (N\Q)] + 2 x k, (5) 

where k is the number of independent parameters and 
ln[P( fc ) (7V|0)] corresponds to log-LF for quantum state 
Pq\ Here is defined as the quantum state that is 
contracted within a k-dimensional parameter space. To 
estimate the two-qubit system, we use Pq 6 \ Pq 5 \ P% 2 \ 
and Pq , which represent the rank-4, 3, 2, and 1 density 
matrices, respectively Among these four states, 

the one that attains the smallest AIC can be regarded as 
the most appropriate one for the following reason. Ap- 
proximation (Q), used for explaining MLE, needs a cor- 
rection depending on the number of parameters || ; that 
is, 

OO OO 

£...]TP o (A)ln[p( fe )(A|0(A))]| t=1 

m= nie— 

»iln[pW(JV|9(JV))]-*. (6) 

Taking this correction into account, the relative entropy, 
Eq. (^|), can be minimized by minimizing the value, 
^rAIC^ k ' (0). Therefore, if the quantum state that min- 
imizes the AIC is chosen from several alternatives, this 
state is the one closest to the true one from the view- 
point of relative entropy (minimum AIC estimate or 
MAICE flj). 

Figure g (b) shows the resulting average Bures dis- 
tances between the true states and the states estimated 
using MAICE. The Bures distances nearly reach the 
lower bounds given by Eq. (|) in the asymptotic regime 
(i.e., in the regime where Ct > 10 3 ), even when estimat- 
ing degenerate states such as APSS and HES (for which 
the AIC mostly gave the rank-2 state as the most suit- 
able one). Therefore, by using MAICE and thus elimi- 
nating the redundant parameters, we demonstrated that 
the CR bound on the statistical errors in quantum-state 
estimation is achievable in the asymptotic regime [ p^[ . 
This means that from give experimental data, we can 
obtain more detailed information about quantum system 
by using MAICE than by any other known estimation 
strategy. 

Finally, we note that the results shown in Fig. |^ (b) 
reveal that in the asymptotic region, the average Bu- 
res distance (i.e., the error) for the entangled state, i.e., 
HES, is the largest of the three states, in spite of the 
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small entropy |IJ] . This may stem from the nature of to- 
mographic measurements using only separable projectors 
for the estimation. Further, we numerically inspected the 
tomographic measurements using inseparable projectors 
and found that it improved the accuracy not only in es- 
timating HES, but also in estimating VNMS JPiJ. The 
latter improvement corresponds to non-locality without 
entanglement |p],p3[. The new strategy, MAICE, should 
make it experimentally possible to analyze such peculiar 
features of quantum mechanics. 

In conclusion, we have demonstrated that by using a 
new estimation strategy based on the AIC, that elimi- 
nates superfluous parameters, the lower bound on sta- 
tistical errors in estimating the quantum states can be 
achieved. This strategy is versatile and can be applied 
to other classes of multi-parameter estimation problems 
for quantum states or even quantum channels. 
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